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High momentum hadron suppression is considered to be an excellent probe of jet- 
medium interactions in QCD matter created in ultra-relativistic heavy ion collisions. 

We previously showed that our dynamical energy loss formalism can accurately ex¬ 
plain suppression measurements at 200 GeV Au-|-Au collisions at RHIC and 2.76 
TeV Pb-|-Pb collisions at the LHC. With the upcoming LHC measurements at no¬ 
tably higher collision energies, there is a question of what differences, with respect to 
the current (2.76 TeV) measurements, can be expected. In this paper we concentrate 
on heavy flavor suppression at the upcoming 5.1 TeV Pb-|-Pb collisions energy at the 
LHC. Naively, one would expect a notably (~ 30%) larger suppression at 5.1 TeV 
collision energy, due to estimated (signihcant) energy loss increase when transition¬ 
ing from 2.76 to 5.1 TeV. Surprisingly, more detailed calculations predict nearly the 
same suppression results at these two energies. We show that this unexpected result 
is due to an interplay of the following two effects, which essentially cancel each other: 
i) flattening of the initial distributions with increasing collision energies, and if) sig¬ 
nificantly slower than naively expected increase in the energy loss. Therefore, the 
obtained nearly the same suppression provides a clear (qualitative and quantitative) 
test of our energy loss formalism. 

I. INTRODUCTION 

High energy heavy flavor suppression [l| is considered to be an excellent probe of QCD 
matter created in ultra-relativistic heavy ion collisions at RHIC and LHC. One of the major 
goals of these experiments is mapping the QGP properties, which requires comparing avail- 
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able suppression data with the theoretical predictions { 2 ^^. Such comparison tests different 
theoretical models and provides an insight into the underlying QGP physics. Having this in 
mind the upcoming 5.1 TeV Pb+Pb measurements at LHC (expected at the end of 2015) 
- and their comparison with theoretical predictions - will provide an additional important 
insight in the jet-medium interactions in QGP created in such collisions. With this moti¬ 
vation in mind, the goal of this paper is providing the heavy flavor suppression predictions, 
and physical interpretation behind the obtained results, for the upcoming high-luminosity 
experimental data at 5.1 TeV Pb-|-Pb collisions at LHC. In particular, we aim assessing 
the differences in the predicted suppression with respect to the already available 2.76 TeV 
measurements at LHC and compare the results of state-of-the-art calculations with simple 
expect at ions/est imat es. 

To generate the theoretical predictions we will use our recently developed dynamical en¬ 
ergy loss formalism, which includes: H dynamical scattering centers, ii) QCD medium of a 


finite size 


, iii) both radiative j^, l6| anff collisional [7[ energy losses, iv) finite magnetic 


mass effects jSj] and v) running coupling j^. This energy loss formalism is based on the 
pQCD calculations in finite size optically thin dynamical QCD medium, and has been incor- 
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porated into a numerical procedure [9| that allows generating state-of-the art suppression 


predictions. The mode 


averaged observables 


has shown to be successful in explaining a wide range of angular 


12| at both RHIC and LHC. Since the angular averaged Raa^ are 


largely insensitive to the medium evolution, angular averaged Raa can be considered an ex¬ 


cellent probe for jet-medium interactions 
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15|; consequently, the suppression predictions 


at 5.1 TeV Pb-|-Pb collisions at LHC, and their comparison with the measurements, will 
allow further testing of our energy loss formalism. 


II. OVERVIEW OF THE COMPUTATIONAL FRAMEWORK 

For generating the suppression predictions, we use the computational procedure from jo]. 
The main features are briefly summarized below, while the full account of the procedure is 
provided in js]. 

The quenched spectra of heavy flavor observables are calculated according to the generic 
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pQCD convolution: 

Efd^a Eid?a{Q) 

dpi dpf 


P{E, ^ Ef) ^ DiQ ^ Hq) ® /{Hq ^ e, J/V^). 


( 1 ) 


In the equation above subscripts ”i” and ”f’ correspond, respectively, to ’’initial” and 
’’final”, and Q denotes heavy quarks. Eid^a{Q)/dpf denotes the initial heavy quark spec¬ 
trum, which is computed at next to leading order according to Q, Hi- P{Ei —)■ Ef) is the 
energy loss probability; this probability includes both radiative and collisional energy 


OSS 


in a finite size dynamical QCD medium, multi-gluon [181 and path-length fiuctuations [1^ 
and running coupling j9|. D{Q ^ Hq) is the fragmentation function of heavy quark Q to 


hadron Hq, where for D and B mesons we use, BCFY 


20| and KLP 2l| fragmentation 


functions, respectively. Finally, decay of B mesons to experimentally measured non-prompt 
J/'ip is represented by /{Hq ^ J/-^) and obtained according to jl6|. 

The expression for the radiative energy loss in a finite size dynamical QCD medium is 
extracted from Eq. (10) in while the collisional energy loss is extracted from Eq. (14) 


m 


0 . 


Path length distributions are taken from 


22 |. 


The angular averaged Raa is a clear jet-medium interaction probe, i.e. it is not sensitive 
on the details of the medium evolution H, 141, so we model the medium by assuming 
constant average temperature of QGP. To determine the average temperatures at 0-10% 
most central collisions, we start from T = 304 MeV (the effective temperature extracted by 


ALICE 


23| for 0-40% centrality), and use the procedure outlined in [l2| (based on gluon 


rapidity density) to determine the temperatures at central collisions at 2.76 and 5.1 TeV 
Pb-I-Pb collisions; for 2.76 TeV 0-10% centrality, this leads to the average temperature of 
313 MeV. To determine the temperature at 5.1 TeV, note that it is expected that the gluon 
rapidity density will be 25% higher at 5.1 TeV than at 2.76 TeV in Pb-|-Pb collisions at 
24l |. Since the temperature is proportional to the gluon rapidity density, i.e. T ~ 


LHC 


{dNg/dy)^P, this leads to ~ 7% higher temperature at 5.1 TeV compared to 2.76 TeV at 
the LHC, i.e. 335 MeV for 0-10% central 5.1 TeV Pb-|-Pb collisions at LHC. Note that, in 
our energy loss calculations, this is the only parameter that differs between the two systems; 
i.e. all the other parameters that enter in the calculations (stated in the next paragraph) 
are the same for the two systems, and correspond to the standard literature values (i.e. no 
parameters are determined through fitting the data). 

The following parameters are used in the numerical calculations: QGP with effective light 
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quark flavors nf = 3 and perturbative QCD scale of Aqcd = 0.2 GeV. The Debye mass is 
taken to be /ig ~ 0.9 (/ig ~ 0.97) GeV for 2.76 (5.1) TeV collision energy, and is obtained by 
self-consistently solving Eq. (7) in [26|. The value for magnetic to electric mass ratio 


is extracted from non-perturbative calculations j27H30| 0.4 < pM/hs < 0.6; the gluon mass 
is Trig = M|, while the charm and the bottom mass are, respectively, M = 1.2GeV 

and M = 4.75 GeV. Path-length distribution, parton production, fragmentation functions 
and decays, which are used in the numerical calculations, are specified above. 


III. RESULTS AND DISCUSSION 


To get an insight of what results we expect at 5.1 TeV collisions at the LHG, we will Erst 
provide a simple analytic estimate for heavy flavor suppression at this collision energy. For 
that purpose, note that radiative energy loss is widely considered to be a dominant energy 
loss mechanism in QGP, so we will use only the radiative contribution for the estimate. 
Since it is also widely assumed that radiative energy loss is proportional to (see e.g. j^l. 
one can estimate that the energy loss at 5.1 TeV should be ~ 25% higher than at 2.76 TeV. 
Based on this, and if we assume that initial distributions can be approximated by power low 
distributions, i.e. da/dp^ ~ l/p”, we can make an estimate on how much larger/smaller 
suppression one would expect at 5.1 TeV compared to the already observed results at 2.76 
TeV. 


It was previously shown that, for radiative energy loss and power low initia 
suppression can be roughly estimated by using the following simple formula 


distributions. 


18 


(1 


lAE. 


(n-2) 


2 E 


(1 


n 


2 AE. 


E 


( 2 ) 


where AE/E is the fractional energy loss. If we assume that, at 2.76 TeV, typical fractional 
energy loss for charm is ~ 30% and for bottom ~ 15%, and that charm and bottom distri¬ 
butions do not notably change between these two collision energies, with n ~ 6.5 (n ~ 6) for 
charm (bottom), the above estimate will straightforwardly lead to the expectation of ~ 30% 
(~ 10%) larger suppression for charm (bottom) at 5.1 TeV compared to 2.76 TeV Pb-|-Pb 
collisions at the LHG. 

Gontrary to these expectations. Figure [U shows that our suppression calculations - ob¬ 
tained from the energy loss formalism outlined in the previous section - provide substantially 
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FIG. 1: Comparison of Raa predictions for heavy flavor at 2.76 and 5.1 TeV. D mesons 
and non-prompt J/'ip suppression predictions, as a function of transverse momentum, are shown 
on the left and the right panel, respectively. Full (dashed) curves correspond to Raa predictions 
at 5.1 TeV (2.76 TeV) Pb+Pb collisions at the LHC. On each panel, the gray bands correspond 


to the finite magnetic mass case (i.e. 0.4 < < 0.6 


27H30|), where the lower and the upper 


boundary correspond, respectively, to hm/ij-e = 0.4 and = 0.6. 


different predictions. From this figure, we actually do not observe any suppression increase 
between 2.76 to 5.1 TeV collisions at the LHC. That is, we obtain the same suppression 
patterns for both charm and bottom probes (D mesons and non-prompt J/'ijj) at these two 
collision energies. This then leads to the question of, why the increase in the collision en¬ 
ergy by almost a factor of 2, leads to the same predicted suppression patterns between the 
two collisional energies, despite the estimated significant (i.e. ~ 30% for charm, see above) 
increase in the suppression? 

To address this question, in Fig. [2]we first compare charm and bottom initial distributions 
between these two collision energies. From this hgure, we see that the distributions at 5.1 
TeV are slightly flatter than at 2.76 TeV, for both charm and bottom, which will have the 
tendency to somewhat lower the suppression at 5.1 TeV compared to 2.76 TeV. Note that 
only the shape of the distributions contributes to the suppression predictions, and from 
Fig. 121 one can observe that the differences in the shape of the distributions are not large. 
Still, this difference in the distributions has a notable (though again not large, i.e. ~ 5%) 
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FIG. 2: Comparison of momentum distributions for charm and bottom at 2.76 and 5.1 

TeV. Charm and bottom momentum distributions, as a function of transverse momentum, are 
shown on the left and the right panel, respectively. On each panel full (dashed) curve corresponds 
to the momentum distribution at 5.1 TeV (2.76 TeV) Pb+Pb collisions at the LHC. 


effect on the suppression predictions, as can be seen in the left panel of Fig. [3l therefore, it 
should be taken into account in the suppression calculations. 

However, what we further see from the right panel of Fig. [3] is that the effect on the 
suppression coming from the energy loss increase between 2.76 and 5.1 TeV (due to the 
increase in average temperature) is also notable but not large, i.e. it corresponds to 5 and 
10%. That is, the energy loss effect on the suppression has about the same magnitude, but 
an opposite direction, compared to the effect of different distributions between these two 
collision energies. The hrst question that we want to address is why the effect of the energy 
loss increase on jet suppression is not larger between these two collision energies, at least 
not for charm quark. That is, based on the common assumption, we have estimated that 
the energy loss increase should be on the order of 25%, which should, therefore, have a more 
prominent (estimated 30%) effect on the suppression. 

Regarding the estimate for the radiative energy loss, note that, while widely used, this 
estimate does not have to be justihed. That is, from Eq. (10) in jsl, which shows the radiative 
energy loss expression in a hnite size dynamical QCD medium, it can be straightforwardly 
observed that the expression nontrivially depends on T. That is, while one can recover a part 
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FIG. 3: Relative increase in Raa between 2.76 and 5.1 TeV. The left panel shows momentum 
dependence of the relative increase in Raa between 2.76 and 5.1 TeV Pb+Pb collisions at the LHC 
due to differences in the distributions; the energy loss is kept fixed and calculated for 2.76 TeV 
case. The right panel shows momentum dependence of the relative increase in Raa between 2.76 
and 5.1 TeV collisions at the LHC due to differences in the energy loss; the momentum distribution 
is kept fixed and calculated for 2.76 TeV case. On each panel, curves that correspond to charm 
and bottom are marked by c and b, respectively, and the magnetic mass is fixed to fJ-M/fJ-E = 0.4. 


with explicit dependence on in this expression, the rest of the expression also depends on 
T, where this extra term considerably modifies the temperatnre dependence. Additionally, 
the collisional energy loss effect, while smaller compared to the radiative, is still important, 
and this effect also has to be taken into account in the suppression calculations. Note 
that, for the collisional energy loss, it is commonly assumed that it has a quadratic (T^) 
dependence on the temperature. However, similarly to the above discussion for the radiative 
energy loss, Eq. (14) from shows a nontrivial temperature dependence, so we will below 
also test whether this simple (T^) assumption is justihed. 

With these goals, in Fig. 01 we plot the relative heavy flavor energy loss increase between 
2.76 and 5.1 collision energies at the LHC. Figures also contain dashed horizontal lines, which 
represent what would be the energy loss increase, if it would indeed have or depen¬ 
dence. For radiative energy loss, we see that, contrary to the common expectations, energy 
loss increase is far from dependence; i.e. it is between linear (for low jet energy regions) 











FIG. 4: Relative energy loss increase between 2.76 and 5.1 TeV. All the panels show 
the momentum dependence of the relative energy loss increase between 5.1 and 2.76 TeV Pb+Pb 
collisions at the LHC. The left, the central and the right panel correspond, respectively, to the 
radiative, collisional and total energy loss case. On each panel curves that correspond to charm 
(bottom) are marked by letter c (b) and the magnetic mass is fixed to hm/ij^e = 0.4. Dashed 
gray horizontal lines represent the energy loss increase, if it would have linear, quadratic or cubic 
temperature dependence. 


and quadratic (which can be reached for asymptotically high jet energies). Consequently, 
for the high momentum heavy flavor hadrons that will be studied at these two collision 
energies at the LHC, the expected energy loss increase is notably smaller than quadratic, 
i.e. it is in the region between 5 — 10% (note that the average temperature increase between 
these two collision energies is ~ 7%). For the collisional energy loss, we also see that energy 
loss increase is far from quadratic, i.e. the increase of ~ 8.5% is constant with momentum 
and it has slightly larger than linear dependence on temperature. Consequently, contrary to 
the common expectation, the total energy loss has also a modest temperature dependence, 
which is close to linear, i.e. between 6 and 10% depending on the jet momentum. This 
modest energy loss increase between these two collisional energies consequently leads to a 
modest increase in the suppression which we observe in the right panel of Fig. [Hi 

Finally, in Fig. O we study the combined effect of the differences in the distributions and 
the energy loss on jet suppression. On the two left panels, we see the effect of the difference 
in the distributions on the jet suppression, while the energy loss is kept hxed. On the two 
central panels, we keep the same distribution, but change the energy loss, while in the two 
right panels both the distributions and the energy loss are changed between the two collision 
energies. From the panels, we see that, while the change in the distribution has the tendency 
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FIG. 5: Analysis of the heavy flavor suppression between 2.76 and 5.1 TeV. The full 
curves correspond to Raas with both the energy losses and the distributions calculated at 5.1 
TeV collision energy. The dashed curves correspond to Raas with both the energy losses and the 
distributions calculated at 2.76 TeV collision energy. The dotted curves correspond to Raas where 
the energy losses are calculated at 2.76 TeV collision energy, while the distributions are calculated 
at 5.1 TeV collision energy. Upper (lower) panels correspond to the charm (bottom) quark. Left 
panels show how the flatter distributions at 5.1 TeV lower the heavy flavor suppression compared 
to the 2.76 TeV case. Central panels show how increase in the energy loss at 5.1 TeV increase the 
suppression compared to the 2.76 TeV case. Right panels show how the above two effects cancel, so 
as to reproduce almost the same suppression at 2.76 and 5.1 TeV Pb+Pb collision energy. On each 
panel, lower (upper) set of curves correspond to the magnetic to electric mass ratio of HM/fJ-E = 0.4 
( ij-m/ij-e = 0 . 6 ). 


to reduce the suppression, the energy loss increase increases the suppression for about the 
same amount, so that the resultant suppression at 5.1 TeV collision energy is almost the 
same as at 2.76 TeV. 

The above obtained numerical result can also be directly estimated from Eq. [2l For 
this purpose, we will take that the energy loss between 2.76 and 5.1 TeV collision energy 
increases by factor 77 , where from Fig. 01 we see that rj ^ 10% for both charm and bottom. 
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Additionally, we will take that the power factor in the initial parton distributions decrease 
by 6] by htting the power low to the ratio of the momentum distributions in Fig. |2l we 
obtain 6 ~ 0.4. By applying these factors into Eq. HI one can straightforwardly obtain 

1 A F! 

Raa{5.1 TeV) ^ i?AA(2.76 TeV) + -^{6 - rj{n - 2)), (3) 

where for 6 and rj estimated above the second additive in the above becomes close to zero. 
Consequently, this estimate also recovers the conclusion of the same heavy flavor suppression 
at 2.76 and 5.1 TeV Pb+Pb collision energies at the LHC. 

IV. CONCLUSION 

In this paper, we provided heavy flavor suppression predictions for the upcoming 5.1 
TeV Pb+Pb collisions at the LHC. Based on our energy loss formalism, we predict the 
same heavy flavor suppression patterns for 2.76 and 5.1 TeV collision energies. This result 
is surprising since, based on the commonly used assumption, a notable increase of the 
suppression is expected at the higher collision energy. We showed that the same suppression 
is a consequence of the interplay between the following two effects: i) a decrease in the 
suppression due to flattening of the initial momentum distributions, and ii) an increase 
in the suppression - though more moderate than expected - due to higher energy loss. 
Consequently, this unexpected, but simple, suppression prediction provides a direct (both 
quantitative and qualitative) test of our understanding of the medium interactions in QCD 
medium created in these collisions. 
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